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Abstract. 

Realistic nucleon-nucleon (NN) interactions, derived within the framework 
of meson theory or more recently in terms of chiral effective field theory, 
yield new possibilities for achieving a unified microscopic description of 
atomic nuclei. Based on spectral distribution methods, a comparison of 
these interactions to a most general Sp(4) dynamically symmetric interac- 
tion, which previously we found to reproduce well that part of the inter- 
action that is responsible for shaping pairing-governed isobaric analog + 
states, can determine the extent to which this significantly simpler model 
Hamiltonian can be used to obtain an approximate, yet very good descrip- 
tion of low-lying nuclear structure. And furthermore, one can apply this 
model in situations that would otherwise be prohibitive because of the size 
of the model space. In addition, we introduce a Sp(4) symmetry breaking 
term by including the quadrupole-quadrupole interaction in the analysis 
and examining the capacity of this extended model interaction to imitate 
realistic interactions. This provides a further step towards gaining a bet- 
ter understanding of the underlying foundation of realistic interactions and 
their ability to reproduce striking features of nuclei such as strong pairing 
correlations or collective rotational motion. 



1 Introduction 

Spectral distribution theory |[T] |2) is an excellent framework for comparing 
the overall behavior of microscopic interactions and uncovering fundamental 
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properties of realistic NN potentials as well as derivative effective interactions 
lIUElE]. Likewise, this information can be propagated beyond the defining two- 
nucleon system to nuclei with larger numbers of nucleons [2| and for higher 
values of isospin [5 1. We search for the level of respect for selected underlying 
symmetries [6, 7| such as the Sp(4) symmetry HIED of iso vector (like-particle 
and proton-neutron, pri) pairing correlations plus an isoscalar pn force and the 
SU(3) symmetry iflOl of collective rotational modes. Such symmetry -respecting 
microscopic model interactions can be used to probe the pairing and rotational 
characteristics of a realistic interaction IfTTl [12l [131 [T4l . which will reflect the 
characteristic properties of the pairing (quadrupole) model Hamiltonian if both 
interactions strongly correlate. As these symmetries are clearly important for 
certain spectral features (for example, observed pairing gaps and enhanced elec- 
tric quadrupole transitions), we have a tool for rapidly assessing the likely suc- 
cess of these interactions for reproducing those spectral features. 

Recent applications of the theory of spectral distributions also include quan- 
tum chaos, nuclear reactions and nuclear astrophysics with studies on nuclear 
level densities, transition strength densities, and parity/time-reversal violation 
(for example, see £15)). The significance of the method is related to the fact 
that low-order energy moments over a certain domain of single-particle states, 
such as the energy centroid of an interaction (its average expectation value) and 
the deviation from that average, yield valuable information about the interaction 
that is of fundamental importance (4j Q~2] [14] [16] \T7] HE) without the need for 
carrying out large-dimensional matrix diagonalization and with little to no lim- 
itations due to the dimensionality of the vector space. Note that if one were to 
include higher-order energy moments, one would gradually obtain more detailed 
results that, in principle, would eventually reproduce those of a conventional mi- 
croscopic calculations. 

We compare three modern interactions, namely, the CD-Bonn [19], CD- 
Bonn+3terms lEOl and GXPF1 f2D . based on realistic nucleon-nucleon po- 
tentials, as well as two model interactions with pairing and quadrupole terms, 
which typically dominate in nuclei. CD-Bonn is a charge-dependent one- 
boson-exchange nucleon-nucleon (NN) potential that is one of the most ac- 
curate in reproducing the available proton-proton and neutron-proton scattering 
data. We employ the two-body matrix elements of the effective interaction de- 
rived from CD-Bonn for Ofrui no-core shell model (NCSM) calculations in the 
fp shell. In addition, the CD-Bonn+3terms interaction introduces phenomeno- 
logical isospin-dependent central terms plus a tensor force with strengths and 
ranges determined in no-core Ohuj shell model calculations to achieve an im- 
proved description of the A = 48 Ca, Sc and Ti isobars. The GXPF1 effective 
interaction is obtained from a realistic G-matrix interaction based on the Bonn- 
C potential [22] by adding empirical corrections determined through systematic 
fitting to experimental energy data in the fp shell. 

Several detailed reviews of the nuclear shell model and its applications have 
been published recently J23] [24] [25] that delve into related key physics issues 
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that we explore. However, the present study is novel and includes /p-shell in- 
teractions, which have been developed since those reviews were completed. 

2 Theoretical Framework 

The theory of spectral distributions (or statistical spectroscopy) is well docu- 
mented in the literature 0] |2] [5] [161 126] |27) . Spectral distribution theory com- 
bines important features, the most significant of which are as follows: 

1 . The theory provides a precise measure, namely, the correlation coefficient, 
for the overall similarity of two interactions. Literally the correlation co- 
efficient is a measure of the extent to which two interactions "look like" 
(are correlated with) one another. In this respect, correlation coefficients 
can be used to extract information how well pairing/rotational features are 
developed in realistic interactions, which may differ substantially from an 
individual comparison of pairing/quadrupole interaction strengths (6). 

2. It gives an exact prescription for identifying the pure zero- (centroid), 
one- and two-body parts of an interaction under a given space partitioning. 
Therefore, major properties follow: 

(a) The correlation coefficients are independent of the interaction cen- 
troids. (A direct comparison of two-body matrix elements provided 
by NN potentials may be misleading, especially when the averages 
of the interactions differ considerably.) 

(b) The pure one-body part of an interaction, the so-called induced 
single-particle energies, is naturally identified in the framework of 
spectral distribution theory and is indeed the average monopole in- 
teraction (compare to ETTO . As such it influences the evolution of 
the shell structure, shell gaps and binding energies [28]. 

(c) The pure two-body part is essential for studies of detailed property- 
defining two-body interactions beyond strong mean-field effects. 

3. The correlation coefficient concept can be propagated straightforwardly 
beyond the defining two-nucleon system to derivative systems with larger 
numbers of nucleons [2 J and higher values of isospin [5 |. This, in addi- 
tion to the two-nucleon information provided by alternative approaches 
(e.g., 1 29 1), yields valuable overall information, without a need for carry- 
ing out extensive shell-model calculations, about the universal properties 
of a two-body interaction in shaping many-particle nuclear systems. 
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For a scalar a = n (isospin-scalar a — n,T) spectral distributiorfj] the cor- 
relation coefficient between two Hamiltonian operators H and H' is defined as 

((fft _ (gt)g)(g/ _ (H') a )) a _ (H^HT - (Hl) a (H') a n . 

CTH&H' OHOH' 

where the "width" of H is the positive square root of the variance, 

(a a H ) 2 = ((H - (H) a ) 2 ) a = (H 2 ) a - ((i/) Q ) 2 , (2) 

and (•••)" denotes an average value related to the trace of an operator divided 
by the dimensionality of the space. The significance of a positive correlation 
coefficient is given by Cohen [30 1 and later revised to the following table: 



Table 1. Interpretation of a correlation coefficient. 

trivial small medium large very large nearly perfect perfect 

0.00-0.09 0.10-0.29 0.30-0.49 0.50-0.69 0.70-0.89 0.90-0.99 1.00 



From a geometrical perspective, in spectral distribution theory every interac- 
tion is associated with a vector and the correlation coefficient £ (Eq. Q]) defines 
the angle (via a normalized scalar product) between two vectors of length a (Eq. 
|2]). Hence, Ch.H' gives the normalized projection of H onto the H' interaction 
(or W onto H). In addition, (C,h,H') 2 gives the percentage of H that reflects the 
characteristic properties of the H' interaction. 

3 Understanding the Nuclear Interaction in Many-nucleon Systems 

We compare the three modern interactions, Hq = {CD-Bonn, CD- 
Bonn+3terms, GXPF1}, based on realistic nucleon-nucleon potentials, and two 
model pairing and quadrupole isoscalar interactions in the fp region by means 
of the theory of spectral distributions. The first model interaction, H S pu), is 
the most general Sp(4)-dynamically symmetric interaction for a system of n va- 
lence nucleons in a 4f2-dimensional space O [9] with two-body antisymmetric 
JT-coupled matrix elements for {r < (s,t); t < u} orbits 0, 

W r J s T tu = -Go^p^jTMoxA-a*. - {-£o[R T + f ] + C}S rt 6 su , (3) 

where G = G + f , E = + g), G, F, E, D and C are interaction 
strength parameters (see Table I in Ref.[9| for parameter estimates). The sp(4) 
algebraic structure is exactly the one needed in nuclear isobaric analog + 
states to describe proton-proton (pp), neutron-neutron (tin) and proton-neutron 

'For n particles distributed over 4Q single-particle states, scalar (or isospin-scalar) distribution 
is called the spectral distribution averaged over the ensemble of all n-particle states (of isospin T) 
associated with the U (4S7) group structure (or U (2f2) (g> U (2)y). 
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(pn) isovector pairing correlations [the first term in (0)], while the second term 
includes isoscalar pn force related to the Eq isospin symmetry term. 

In addition, we construct an extended pairing plus quadrupole model, Hm, 
by including an additional traceless quadrupole-quadrupole interaction that be- 
ing symmetric under SU(3) breaks the Sp(4) symmetry, 

Hm = H sm + H^{2), H Q = -|q ■ Q, (4) 

where \ is the only parameter in the present analysis and is determined by opti- 
mum correlation coefficients. The Hq (2)-term is the part of the pure two-body 
Hq(2) interaction that is not contained in (is orthogonal to) the Sp(4) interac- 
tionj. This is because the Sp(4) interaction contains a nonnegligible part of the 
quadrupole-quadrupole interaction that is revealed by the correlation between 
Hq and H Bp ^y Namely, in the scalar case it is 15% (I/7/2), 29% (I/5/2) and 
29% (2pi/22p3/2), and for the T=l part of the interactions in the isospin-scalar 
case, it is 34% (l/ 7/2 ), 58% (l/ 5/2 ) and 58% (2p 1/2 2p 3/2 ). This is probably 
one of the reasons why the Sp(4) model interaction turns out to work rather well 
despite no explicit appearance of the quadrupole-quadrupole interaction. 

The present study focuses on the weaker but property-defining two-body 
part of the interactions, such as H (2), in the I/7/2 @ and upper ^p- shell Q 
domains. Such a partitioning of the fp oscillator shell follows naturally from a 
splitting of these two regions by a strong spin-orbit interaction. Several corre- 
lation coefficients are of particular interest. Specifically, the overall similarity 
between a realistic interaction and the extended pairing +quadrupole model in- 
teraction in the (isospin-) scalar case is estimated by the ChJ(2) h m corre lation, 
while the capability of a realistic interaction to describe rotational collective mo- 
tion, and hence to reproduce rotational bands and enhanced electric quadrupole 
transitions, can be detected via its correlation to the full Hq(2) quadrupole- 
quadrupole two-body interaction, ChJ(2) h q (2)- isospin-scalar space par- 
titioning is where the ability of a realistic interaction to form correlated pairs 
and hence reproduce prominent pairing gaps is detected via Cjj D (2) h 4 • ^ n tne 

present analysis, the ChJ(2) h q (2) anc * ^h^(2) h (4) correlation coefficients are 
independent of the quadrupole/pairing interaction strength parameters. 

In the detailed case of isospin-scalar distribution in the I/7/2 orbit, the CD- 
Bonn, CD-Bonn+3terms and GXPF1 interactions are found to contain on av- 
erage 59%, 77%, and 78%, respectively, of the pairing+quadmpole interaction. 
This percentage goes up to 91%, 97%, and 92%, respectively, for the highest 
possible isospin group of states for all the nuclei with valence protons and neu- 
trons occupying the I/7/2 shell (Fig. [T]). For these states, the strongest correla- 
tion was observed between the CD-Bonn+3terms and the pairing+quadrupole 



'Such a Hamiltonian (Eq.[4) does not affect the centroid of -ff s p(4) because Hq(2) is traceless 
and hence preserves the shell structure that is built into ^ s p(4) a nd established by a harmonic 
oscillator potential and as a result is favored in many studies [12 14. 31|. 
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model interaction, where other types of interactions accounted for in the realistic 
interaction represent only 3% of it. They constitute 8% of the GXPF1 inter- 
action, and 9% of CD-Bonn. While both interactions, CD-Bonn+3terms and 
GXPF1, exhibit a well-developed pairing character compared to CD-Bonn, the 
latter appears to build up more (less) rotational collective features that are out- 
side of the scope of the T = 1 (T = 0) Sp(4) interaction. 




Figure 1. Geometrical representation of the T = 1 CD-Bonn (light blue), CD- 
Bonn+3terms (red) and GXPF1 (green) interactions, in an abstract operator space, where 
the horizontal plane is spanned by the orthogonal linear operators, the pure two-body 
H s p(4) and Hq(2) model Hamiltonians, both linearly independent of the residual inter- 
action operators represented by the vertical axis. The orientation of the vectors remains 
the same for any particle number n > 2 and for all T — n/2 cases. 

The abovementioned strong correlation coefficients, H J , imply that 
both realistic and Hm interactions are expected to yield energy spectra of a sim- 
ilar pattern. Indeed, for these cases the pairing +quadrupole model interaction 
appears to be a very good approximation that provides a reasonable description 
of the energy spectra of the nuclei in the l/ 7 / 2 level (Fig. |2j. 



CD-Bonn CD-Bonn+ GXPF1 
3terms 

>-<E -£•! 



42 Ca 4 -Sc 42 Ti 54 Co 54 Fe 



Figure 2. Energy spectra of T = 1 states predicted by the CD-Bonn (light blue), CD- 
Bonn+3terms (red) and GXPF1 (green) interactions. Each is compared to the model 
Hamiltonian Hm (black) with \ — 0.071, 0.036 and 0.055, respectively. For compari- 
son, the experimental T — 1 energy spectra of the A — 42 Ca, Sc, Ti isobars (blue) and 
A = 54 Co and Fe isobars (magenta) are also shown. 

In the upper fp domain the outcome yields strong correlation of the pure 
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two-body GXPF1 interaction with the pairing+quadrupole extended model in 
both scalar and isospin-scalar distributions (Fig. [3]). Particularly, the outcomes 



(CD-Bonn. Sp(4n 
(CD-Bonn, H ^(2)) 
(CD-Bonn, Hi 



yywwwwu 



AAATVA 




Figure 3. Correlation coefficients of the pure two-body CD-Bonn (blue), CD- 
Bonn+3terms (red) and GXPF1 (green) interactions with H sv u\ (filled symbols), 
Hq{2) (transparent symbols) and Hm (empty symbols) in the upper fp shell for the 
isospin-scalar distribution. For each valence particle number, n, the isospin T varies as 
j, ? — 1, . . . , 0(|). The figures are symmetric with respect to the sign of n — 2Q. 

show very good scalar-distribution correlations of GXPF1 with the Sp(4) dy- 
namically symmetric interaction and with the Hq (2) quadrupole-quadrupole in- 
teraction (Table|2]i. Rotational features within many-nucleon systems in the up- 
per^? domain are found to be more fully developed for GXPF1 and less for 
CD-Bonn+3terms and CD-Bonn. The isospin-scalar space partitioning demon- 
strates a tendency in GXPF1 towards the formation of correlated pairs in the 
highest possible isospin groups of states (Fig. |3). 

The different extent to which the GXPF1 interaction compared to the CD- 
Bonn and CD-Bonn+3terms interactions reflects development of pairing cor- 
relations and collective rotational modes in the upper fp domain may be the 
reason why their pure two-body part do not correlate as strongly as, for exam- 
ple, CD-Bonn and CD-Bonn+3terms do. Namely, in the scalar case the pure 
two-body correlations are 0.90 (between CD-Bonn and CD-Bonn+3terms) and 
only 0.56 (CD-Bonn and GXPF1) and 0.53 (CD-Bonn+3terms and GXPF1). 
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Table 2. Correlation coefficients for many-nucleon systems of the Ho(2) pure two- 
body part of the CD-Bonn, CD-Bonn+3terms and GXPF1 interactions with H S pu) and 
Hq(2), with the pure two-body full quadrupole-quadrupole interaction, Hq(2), and with 
the extended pairing+quadrupole model interaction Hm (Eq.0. 





CD-Bonn 


CD-Bonn+3terms 


GXPF1 


CHo(2),H sp(4) 


0.55 


0.50 


0.65 


Ch (2),H^(2) 


0.14 


0.20 


0.51 


Ch (2),Hq(2) 


0.28 


0.33 


0.67 


Ch (2),H m 


0.57 


0.54 


0.83 



In the isospin-scalar case, the correlations vary slightly with the particle num- 
ber and isospin and they are on average, 0.88 (between CD-Bonn and CD- 
Bonn+3terms), 0.40 (CD-Bonn and GXPF1), and 0.37 (CD-Bonn+3terms and 
GXPF1). In addition, one can compare the significant monopole influence of 
the three interactions, which is very similar for all when averaged over the 
isospin values. However, in the isospin-scalar distribution, the correlation co- 
efficients involving the induced effective one-body contribution differ between 
GXPF1 and the two CD-Bonn interactions. Their behavior, especially below 
mid-shell, reflects the similarity of the corresponding T = induced single- 
particle energies and the opposite signs of the corresponding Aj^ 1 (for 2p 3 / 2 ) 
and AjT^ 1 (for I/5/2 ) pure one-body interactions. 

Individual orbit analysis, including the I/7/2. 1/5/21 2pi/2> an d 2p 3 / 2 lev- 
els, shows considerably stronger correlation of all the interactions with the pair- 
ing+quadrupole model interaction (up to . 8 — 1 . 00) as well as in nuclear systems 
with more than two nucleons. However, more prominent differences among the 
interactions appear in the multi-j upper fp domain especially concerning both 
CD-Bonn interactions. This may indicate that the inter-orbit interactions do not 
respect strongly the symmetries imposed in the model interactions. In addition, 
the interaction strengths may differ from one orbit to another. While they do not 
affect correlation coefficients in the singe- j cases, their relative strength is of a 
great importance for multi-j analysis. In addition, the difference in the behavior 
of CD-Bonn+3terms within both regions, I/7/2 and upper fp, may reflect the 
fact that this interaction was determined through a reproduction of the energy 
spectrum and binding energy of A = 48 I/7/2 nuclei. 

In the upper fp region, the extended Hm pairing+quadrupole interaction 
is strongly correlated with the pure two-body GXPF1 interaction especially 
in the scalar distribution (Table [2]) and for this reason can be used as a good 
approximation. This is reflected in the quite good agreement between the 
experimental low-lying energy spectra of 58 Ni and 58 Cu and the theoretical 
prediction based on the Hm model interaction with x = 0.027 (see Eq. @]l in 
the l/ 5 /22pi/ 2 2p 3 /2l59/2 major shell (Fig. Hi. 
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Figure 4. Theoretical ('th') low-lying energy spectra for 58 Ni (left, blue) and 58 Cu (right, 
red) compared to experiment ('exp', black). The theoretical calculations are performed in 
the l/5/22pi/22p3/2l<?9/2 major shell with the Hm model interaction with \ — 0.027 
and with single-particle energies derived from 57 Ni experimental energy levels. 



In summary, the present study reveals important information, within the 
framework of spectral distribution theory, about the types of forces that dom- 
inate the /p-shell CD-Bonn, CD-Bonn+3terms and GXPF1 interactions in nu- 
clei and their ability to account for many-particle effects such as the formation of 
correlated nucleon pairs and enhanced quadrupole collective modes. The results 
also illustrate that interactions, which were found to strongly correlate, produced 
energy spectra of a similar pattern. 
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